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Abstract. A set of exact closed-form Bloch-state solutions to the stationary Gross- 
Pitaevskii equation are obtained for a Bose-Einstein condensate in a one-dimensional 
periodic array of quantum wells, i.e. a square-well periodic potential. We use these 
exact solutions to comprehensively study the Bloch band, the compressibility, effective 
mass and the speed of sound as functions of both the potential depth and interatomic 
interaction. According to our study, a periodic array of quantum wells is more 
analytically tractable than the sinusoidal potential and allows an easier experimental 
realization than the Kronig-Penney potential, therefore providing a useful theoretical 
model for understanding Bose-Einstein condensates in a periodic potential. 
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1. Introduction 

Bose-Einstein condensates (BECs) in periodic potentials have attracted great interest 
both experimentally and theoretically during the past few years jT], [2] . A major reason is 
that they usually exhibit phenomena typical of solid state physics, such as the formation 
of energy bands 0, 0], Bloch oscillations [5j E], Landau-Zener tunneling [TJ, El El ED] 
between Bloch bands and Josephson effects [HI [12], etc. The advantage of BECs 
in periodic potentials over a solid-state system is that the potential geometry and 
interatomic interactions are highly controllable. Such a BEC system can therefore 
serve as a quantum simulator [13] to test fundamental concepts. For instance, the 
Bose-Hubbard model is almost perfectly realized in the BEC field, hence enabling an 
experimental study of the quantum phase transition between a superfluid and Mott 
insulator [T31IT5]. 

Research so far has been primarily focused on BECs in two types of periodic 
potentials. The first type is the sinusoidal optical lattice [U [2] . Experimentally created 
by two counter-propagating laser beams, the sinusoidal optical lattice consists of only 
one single Fourier component. Most studies on BECs in this type of potential ask for 
the help of numerical simulations since analytical solutions are lacking. By contrast, the 
second one is the so-called Kronig- Penney potential [T7], [TU [19]. In the BEC field, the 
Kronig- Penney potential as shown in Refs. [TFJ [TBI [19] is usually referred to as a periodic 
delta function potential. However, in the original work [20] and the field of condensed 
matter physics [21 J, the Kronig- Penney potential is also used as the periodic rectangular 
potential. To avoid the confusion, we adopt the notion in the BEC field and refer to the 
Kronig- Penny potential as a periodic delta function potential. With this understanding, 
the periodic rectangular potential in this paper is called as a periodic array of quantum 
wells. The Kronig-Penny potential admits an exact solution in closed analytical form, 
leading to general expressions that can simultaneously describe all parameter regimes. 
Nevertheless, it's very difficult to realize a Kronig-Penney potential in experiments. It is 
therefore instructive to seek a periodic potential that not only permits an exact solution 
in closed analytical form, but also is hopeful to be realized experimentally. The search 
for such type of potential is justified by the fact that the fundamental properties of a 
BEC in a periodic potential should not depend on the potential shape [18]. So theorists 
are actually at liberty to select the form of a periodic potential for the convenience of 
their study. 

One such option is provided by a periodic array of quantum wells separated by 
barriers [2~T] . On the experimental side, this potential can be generated by interference 
of several laser beams. Since two interference counter-propagating laser beams form a 
sinusoidal potential that contains one single Fourier component, we expect more Fourier 
components to be involved by using several counter-propagating laser beams. When 
frequencies of these beams are multiples of the fundamental, interference of them would 
result in a periodic array of quantum wells. An experimental scheme to create such 
unconventional optical lattices has recently been demonstrated in Ref. [22]. On the 
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theoretical side, it will be shown in this paper that exact closed-form solutions exist for 
a periodic array of quantum wells. In fact, such potential virtually becomes a Kronig- 
Penney potential, i.e., a lattice of delta functions, in the limit when the width of barriers 
becomes much smaller than the lattice period. We are therefore motivated to launch a 
systematic study on a BEC in a periodic array of quantum wells. 

In this article, we derive a set of exact Bloch-state solutions to the stationary Gross- 
Pitaevskii equation (GPE) for a BEC in a one-dimensional periodic array of quantum 
wells. All our exact solutions, in the limit of varnishing interatomic interaction, are 
reduced to their counterparts in the linear case, i.e. the Bloch states of the stationary 
Schrodinger equation with a one-dimensional periodic array of quantum wells. We apply 
these solutions to analyze the structure of Bloch bands, the compressibility, effective 
mass and the speed of sound as functions of both potential depth and the strength of 
interatomic interaction. Special emphasis is given to the behavior of the compressibility 
and effective mass. 

The outline of the paper is as follows. In Sec. 2, we introduce notations and describe 
basic theoretical framework of our study. In Sec. 3, the general solutions of GPE in a 
single quantum well are derived in detail. In Sec. 4, we investigate the lowest Bloch band 
for a BEC in a periodic array of quantum wells. A comprehensive analysis is presented 
in Sec. 5 that explains the dependence of Bloch band, the compressibility, effective 
mass and the speed of sound on the potential depth and the strength of interatomic 
interaction. Finally, we discuss their experimental implications followed by a summary 
in Sec. 6. 



2. Mean-field theory of Bose-Einstein condensates 

We consider a BEC is tightly confined along the radial directions and subjected to a 
periodic potential in x-direction. The periodic potential V pot (x) is assumed to be a 
periodic array of quantum wells in the form 

+oo 

V pot (x)= V q (x-nT), (1) 

n=— oo 

with 

V q (x) = {\ 0<X ~ a > u (2) 

q K ' | sE R a<x<a + b, w 

where a is the well width and b is the barrier width. In Eq. (CQ), the V po t (x) has a 
periodicity of T = a + b. The s in Eq. ([2]) is a dimensionless parameter that denotes 
the strength of the V pot (x) in units of the recoil energy Er = h 2 q B /2m, with q B = n/T 
being the Bragg momentum. 

We restrict ourselves to the case where the BEC system can be well described by the 
mean-field theory. The parameter characterizing the role of interactions in the system 
is g3D~n, where g 3 n = A.Tih 2 a s /m is the two-body coupling constant and n is the 3D 
average density. Here a s > is the 3D s-wave scattering length. At the mean-field level, 
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descriptions of a BEC system are given by the stationary GPE (or nonlinear Schrodinger 
equation). In our case, the confinement along the radial direction is so tight that the 
dynamics of the atoms in the radial direction is essentially frozen to the ground state of 
the corresponding magnetic trap. As shown in Ref. [23], the effective coupling constant 
can be deduced as g = g 3 o/2ll with l being the length scale of the magnetic trap. In 
this limit, the stationary 3D GPE therefore reduces to a ID equation that reads [2] 

h 2 d 2 



+ V pot (x) +gnT\ip\ ip (x) = pup {x 



(3) 



2m dx 2 

where m is the atomic mass, \x is the chemical potential and the order parameter ip{x) 
is normalized according to J Q T dx\ip(x)\ 2 = 1. 

Despite non-linearity, Eq. ([3$ permits solutions in the form of Bloch waves [3[ [25] 

Mx)=e ikx Mx), ( 4 ) 

where k is the Bloch wave vector and <p k (x) is a periodic function with the same 
periodicity as the V pot (x). We point out that Eq. (T4|) does not exhaust all possible 
stationary solutions of GP Eq. ([3]) due to the presence of the nonlinear term. Except 
for the Bloch-form solutions, the GP Eq. ([3]) with a periodic potential also allows other 
kinds of solutions, for example, period-doubled state solutions [IB] . 

The GP equation ([3]), in terms of the function 4>k{x), can be rewritten as: 

-ihd x + k) 2 



2m 



+ V pot + gnT\(j> k \ 



<p k = n(k)4> 



k ■ 



Note that the chemical potential pi(k), which is derived from Eq. §5§ as 



/*(*) 



dx(p* k 



-ihd x + kf 
2m 



+ V po t + gnT\<p k 



usually does not coincide with the energy e(k) of the BEC system defined by 



e(k) 



dx<f>t 



-ihd x + k)' 



+ Vpot + -gnT\cf) k 



(5) 



(6) 



(7) 



2m puu 2 

Comparison of Eqs. (jSJ) and (J7J) indicates that e{k) equals fi(k) only when interactions 
are absent. Generally, e(k) and fi(k) are related to each other through following 
definition EH [25] 



d [ne(k)] 
dn 



(8) 

by assuming following ansatz [TF], [181 [12] 



Now we seek exact solutions of Eq. 
for the wave function ip(x) 

V> (x) = a/ p (x) exp [-iQ (x)) , (9) 

where the density function p (x) is nonnegative and the phase function (x) is real. 
Substituting Eq. Qj into Eq. ([3]) and re-scaling equations, we obtain 

= 2r/p 3 + 4 fa - Vpot) P 2 - I3p - 4a 2 , (10) 
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and 

= I dx-^—, (11) 

J pyx) 

where the length scale is T/ir, the t] = gnT/E^ represents the nonlinear interaction, 
and /3 and a are integral constants. 



3. General solution in a single quantum well 

We then proceed to solve the GP Eq. (fTUl) in a single quantum well V S i g defined by 



V Sl g (X) 



< x < a, 
a < x < a + b. 



(12) 



As is well known, general solutions for Eq. ([TO]) with a constant potential can be 
expressed in terms of the Jacobi Elliptic functions [2J ETJ EE]- In our case, we derive 
exact solutions to Eq. (flOl) separately in the two regions shown in Eq. (fl2|) . 

In the region < x < a, the V s i g (x) is zero. Hence the exact solutions of Eq. (TTDT) 
have the following general form 

2a 2 



Pi (x) = A 



A 



where 



.4 



A(2K 2 + Arj)_ 



a 



SN 2 {Kx + 6, n\) 



(13) 



7>7 







2K 2 1 AK 2 (2K 2 + Ar]) n 

2 {2AK 2 + A 2 T]) 2 + 8{A V + K 2 ) a 2 



K 2 + A + 



2AK 2 + A 2 r] 

„2 



(14) 



a 



A (2K 2 + An) 

In Eq. ( TTB1 . SN is the Jacobi elliptic sine function and n\ denotes the modulus whose 
range is restricted within [0, 1]. In this general solution, the free variables are the 
translational scaling K, the translational offset 5, and the density offset A. In the limit 
of i] = 0, the solution ffTB"]) can be reduced to 



Pi (x) = A 



A 



a 2 \ 
AK 2 ) 



sin 2 (Kx + S) 



(15) 



with n\ = 0, (3 



8A 2 K 2 +8a 2 
2A 



and /i = K 2 . 



In the region a < x < a + 6, the V S i g (x) is a constant. In this region, Eq. (flOl) 
admits two kinds of exact solutions, depending on whether there is a node within the 
barrier. 

The first type of solutions contains no node within the barrier and has the form: 

2a 2 



p 2 (x) = B 



B 



B (2Q 2 - Br,) 



SC 2 (Qx + 7 ,n 2 ) 



(16) 
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with 



n 



l> 



1 - 



B 



a 



2Q 2 BQ 2 (2Q 2 - Br,) 

2 



s-Q 2 



B 



a 



(3 = ABQ 2 - 2B 2 r) + 



B (2Q 2 - Br]) 
(8Br] - 8Q 2 ) a 2 



(17) 



2BQ 2 - B 2 r] 
In the limit of rj = 0, this solution is reduced to 

„2 



p 2 {x) = B+[B 



a 



BQ< 



sinh (Qx + 7) 



[11 



with n 2 = 1, fi = s - Q 2 and (3 = (AB 2 Q 2 - 4a 2 )/ B. 

The second type of solutions admits only one node within the barrier and is 
expressed as 



B 



with 



P'2 



Tin 



l> 







8Q 2 



Q 2 B Q^BQ 2 - 16a 2 r] 



Br] 



NC 2 (Qx + 1 ,n 2 ), (19) 



1 _ _B_ ^BQ 2 - 16 

2 ~ TOO 1 " 1 + 2VBQ 



a 2 r] 



B Q^/BQ 2 - lQa 2 r] 



i6g 4 ■ 

32Q 2 a 2 \[B\jBQ 2 - 16a 2 r] 



(20) 



B 2Q 
which again can be reduced in the limit of 77 = to: 



P2 0) 



B 



B 



cosh 2 (Qx + 7) 



(21) 



8Q 2 8Q 2 

with n 2 = 1, fx = s - Q 2 and (3 = (B 2 - UQ 2 a 2 )/2B. In Eqs. (JTBJ) and(fT9l). and 
iVC are also the Jacobi elliptic functions with modulus n| and the free variables are B, 
Q and 7. 

Note that all above solutions, in the limit of rj = 0, are nothing but the stationary 
solutions for the linear Schrodinger equation. 



4. Bloch bands and group velocity 

So far we have ignored Bloch wave condition (j3J) and solved the GP equation for a single 
quantum well for specific regions. Next we seek the global solution to the GP equation 
defined on the whole x axis that satisfies the Bloch wave condition. 

Assume that the solutions given in Eqs. (|T3|) . (fTEl) and (jT9l) respectively comprise 
a segment of the complete Bloch wave stationary solution of the GP equation with the 
potential given by Eq. ([1]). We then extend the wave function ipk{x) originally defined on 
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(0, a + b] to the whole x axis and construct the ultimate Bloch wave solution according 
to the Bloch condition: 

Mx + T) = e lkT ^ k {x), (22) 

where k is the Bloch wave vector defined by 

r T dr 

k = -S(x) = -a / (23) 
Jo P( x ) 

Imposing the boundary condition that the ipk{ x ) is continuous at x = a and using Eq. 
( I2"2"j) . we obtain 

Pl (a) = p 2 (a) , (24) 
d x pi (a) = d x p 2 (a) , (25) 
Pl (0) = p 2 (T) , (26) 
d xPl (0) = d xP2 (T) . (27) 

Two additional constraints are the continuity of the chemical potential on the boundary 
and the normalization condition for the wave function, i.e. 

Mi = M2, (28) 

and 

T pa pT 

\?pk(x)\ 2 dx = / P i (x) dx + / P 2 (x) dx — 1. (29) 

Jo J a 

In principle, Eqs. (I23H29P provide a complete set of equations for us to determine the 
unknown parameters A, B, K, Q, 5, 7 and a. Once these parameters are found, we 
obtain the Bloch band of the periodic system. Note that in the region a < x < a + b, 
the first type solution P 2(x) is adopted when k ^ ±qB , while the second type P 2(x) is 
used when k = ±qs- 

In Fig. [TJ the lowest Bloch band and the corresponding group velocity as a function 
of quasi momentum k are presented for various potential depth s and interatomic 
interactions gn. Quantum well width a = 0.6, and barrier width b = 0.4 are used 
throughout this paper. Here, the group velocity is defined by Vk = de/dk. The states 
with k = and k = ±q B in Fig. [1] respectively correspond to the stationary condensates 
at the bottom and top of the lowest Bloch band. The state with k ^ and k 7^ ±qs 
in Fig. [U on the other hand, describe a condensate where all atoms occupy the same 
single-particle wave function and move together in the periodic potential with a constant 
current riVk- 

Furthermore, Figs. [U(a), (e) and (b) demonstrate that when the potential depth s 
is fixed, the interatomic interactions affect the group velocity more conspicuously than 
the Bloch band. Yet for given interatomic interaction gn, Fig. [1] (c) and (d) show 
that the Bloch band becomes more and more flatter with increasing potential depth. 
Eventually, when the potential wells are sufficiently deep, the condensate becomes so 
localized in each quantum well that an adequate description can be obtained by directly 
using the tight-binding model [21 
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Figure 1. (a) and (b) the lowest Bloch band and group velocity with s = 5 for gn = 0, 
0.1-Er and 0.5-Er; (c) and (d) the lowest Bloch band, group velocity with s — 5, 10, 15 
when gn = O.IEr. (e) in order to compare changes of the energy band, we vertically 
translate the Bloch band of gn = Q.IEr and 0.5-Er to the initial position oi gn = 0. 



5. Compressibility, Effective mass and sound speed 

Now we apply our exact solutions to study the compressibility, the effective mass and 
the sound speed of a BEC in a periodic array of quantum wells. 

We start by calculating the compressibility k. In thermodynamics, k is defined as 
the relative volume change of a fluid or solid with respect to a pressure (or mean stress) 
variation. In our case, the compressibility k is given by|24[ [25] 



= n d~n (30) 
For a BEC system with repulsive interatomic interaction, the periodic potential traps 
atoms and enhances the repulsion. A reduced compressibility k is therefore expected. 
We illustrate this point in detail in the following. 

In the uniform case of s = 0, the chemical potential is linearly dependent on the 
density expressed by \i = gn. Thus k" 1 = gn is proportional to the density. When 
s 7^ 0, we substitute the general expression of fi given by Eq. (jSJ) into Eq. ( [30]), and 
obtain the ft -1 for a BEC system in a periodic array of quantum wells. The calculated 
k _1 is plotted in Fig. [2] as a function of the interatomic interaction gn for different 
s. The figure demonstrates that the n~ l increases with s, typical of a wave function 
localized at the bottom of each quantum wells. Compared to the uniform case, the k^ 1 
increases linearly only for small gnjE^. Whereas for large gn/En, the growth of k^ 1 
develops a nonlinear dependence on gn/En. 

We now give an analytical explanation to the behavior of hT 1 shown in Fig. El 
Assume that k~ x is related to s by the following expression [23] when gn/E^ is small 

z^ 1 = 9(s)n, (31) 
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where 



/i = /V=o + 9 (s) n, 



(32) 



in which p gn =o depends on the potential depth, but not on density. The quantity 
g(s) in Eq. (13T]) acts as an effective coupling constant. In case Eq. (}3~Tj) is valid, 
the compressibility of a BEC in a periodic array of quantum wells with g is virtually 
transformed to the compressibility of a uniform BEC with the g{s). Thus by simply 
replacing g by g(s), we can view our system as if there is no periodic potential [21] as 
far as the compressibility is concerned. 

To obtain the form of g{s), we substitute Eq. into Eq. (130]) yielding [24"] 



where 0^=0 is the ground state solution of Eq ([5]) for r\ = 0. Comparison of Eq. (1331) 
with Eq. ( 13TI) gives 



where pi (x) and p2 (x) are solutions respectively in well and barrier. 

We plot the k~ x / gn as a function of the potential depth s for gn = O.IEr and 
gn = 0.5Er in Fig. [3l To compare with the behavior of the effective coupling constant 
g defined by Eq. f l3~l) . the function of g/g with s is also plotted. Fig. [3] shows that the 
linear dependence of n" 1 on gn breaks down. However, with the increasing of s, the law 
of K _1 = g~n becomes to be applicable. 




(33) 




(34) 



which in our formulation has the following form: 




(35) 



2.0 



1.5 - 



s=0 

s=5 

s=10 




1.0 
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Figure 3. n~ 1 /gn for gn = O.IEr (dashed-dotted line) and gn = 0.5Er (short dashed 
line) as a function of the potential depth s, comparing with the effective coupling 
constant g/g (solid line). 



We then consider the effective mass. A BEC trapped in a periodic potential can be 
approximately described by a uniform gas of atoms each having an effective mass m* 
defined by H [25] : 

1 d 2 e(k) 



m* dk 2 



(36) 

k=0 

The dependence of effective mass m* (k = 0) on the potential depth s for gn = 0, 
gn = O.IEr and gn = 0.5Er is demonstrated in Fig. HI According to Fig. [31 when 
s = 0, the effective mass m* is reduced to the bare mass m. Whereas when s increases, 
for example, to s = 30, the m* becomes two orders larger in magnitude than the m. This 
increase of m* with s can be explained by the slow-down of the particles during their 
tunneling through the barriers. Fig. H] also demonstrates that the the m* effectively 
decreases with increasing interactions. This is because repulsion, contrary to the lattice 
potential that serves as a trap, tends to increase the width of the wave function which 
favors tunneling. This is the so-called screening effect of the nonlinearity [5]. 

As is emphasized in Ref. [23], the m* is determined by the tunneling properties 
of the system, thereby exponentially sensible to the behavior of wave function within 
the barriers. Thus any small change in the wave function will significantly affect the 
value of m*. As a result, the conventional Gaussian approximation [23] in the tight- 
binding limit can not be employed to calculate the m*. In this aspect, a periodic array of 
quantum wells as a solvable model, provides a better choice than the sinusoidal potential 
in studying the m* of a BEC in a periodic potential. 

Finally, we proceed to study the sound speed. Sound is a propagation of small 
density fluctuations inside a system [231 [23 [291 USB I5T] . The key point in studying 
sound is to find the sound speed. The speed of sound is important for two simple 
reasons: (i) it is a basic physical parameter that tells how fast the sound propagates 
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in the system, and (ii) it is intimately related to superfluidity according to Landau's 
theory of superfluid. Because of these, the sound propagation and its speed were one 
of the first things that have been studied by experimentalists on a BEC since its first 
realization in 1995 [32, 33J. 

The first step to derive sound speed in a BEC is to find the ground state since 
it acts as a media for the sound propagation. Next, one determines the sound speed 
by perturbing the ground state. Traditionally, there are two equivalent definitions for 
the sound speed [25]. In the first definition, sound is regarded as a long wavelength 
response of a system to the perturbations. Sound speed can be extracted from the lowest 
Bogoliubov excitation energy, which is characterized by the linear phonon dispersion 
with a finite slope. We emphasize that the physical meaning underlying the Bogoliubov 
spectrum is very different form that of the Bloch bands discussed in Fig. 1. The 
Bloch bands refer to states which involve a motion of the whole condensate through the 
periodic potential. However, the Bogoliubov spectrum describes small perturbations 
which involve only a small portion of atoms. The non-perturbed condensate acts as 
a medium through which the perturbed portion is moving. In other words, the Bloch 
band gives the energy per particle of the current states. Being multiplied by N, the 
Bloch band energies obviously exceed the energies of the Bogoliubov excitations. In the 
second definition, the BEC system is viewed as a hydrodynamical system. Accordingly, 
the sound speed in a BEC assumes following standard expression [2U [251 ESS f3Tj 

V sound = ; (37) 

V Km* 

Here we adopt the second definition of sound speed in Eq. (j3Tj) in following calculations, 
using our previous derivation of the compressibility and effective mass. 

The calculation of sound speed as a function of the potential depth s is plotted 
in Fig. [5] for gn = O.IEr and gn = 0.5Er. The figure demonstrates that sound 
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velocity decreases when the potential depth is increased. This can be explained by the 
competition between the slowly decreasing k and the increasing m* when lattice depth 
is increased. 

6. Conclusion 

In typical experiments to date, the relevant parameters are usually chosen as follows: the 
interatomic interaction gn ranges from 0.02Er to IEr [TJ [2]; the depth of the periodic 
potential s can be adjusted from OEr to 12Eji [15J, whereas the BEC system is still 
kept in the superfluid state. In particular, for a one-dimensional periodic potential, the 
transition to the insulator phase is expected to happen for very deep lattice. Thus there 
is a broad range of potential depths where the gas can be described as a fully coherent 
system within the framework of the mean field GPE. Hence the range of parameters 
in our model fit well in the current experimental conditions. Further more, a periodic 
array of quantum wells could be experimentally generated by the interference of serval 
two-counter-propagating laser beams [22J. However, we would like to point out that 
our study is based on GPE. In this mean-field theory, all quantum fluctuations and 
temperature effect are ignored. Thus in order to study the effects of temperature or 
fluctuations, one has to use other theories [3lj, especially near the transition point of 
superfluid and Mott insulator. 

In this paper, we obtain a set of exact closed-form Bloch-state solutions to the 
stationary GPE for a BEC in a one-dimensional periodic array of quantum wells. These 
solutions are applied to calculate the Bloch band, the compressibility, effective mass and 
speed of sound as functions of the potential depth and the interatomic interaction. As a 
result, this type of periodic potential provides a useful model for further understanding 
of BECs. 
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